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ABSTRACT 

The leading order finite-width corrections to the equation of motion de- 
scribing the motion of a domain wall are derived. The regime in which this 
equation of motion is invalid is discussed. Spherically and cylindrically sym- 
metric solutions to this equation of motion are found. We also clarify a miscon- 
ception that has arisen in recent years regarding the rigidity (or otherwise) of 

cosmic strings. 
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There has recently been renewed interest in domain walls as possible sources of large 
scale structure in the universe 1 . With a few exceptions, most of the work done on the 
motion of walls has used the zero thickness approximation. The numerical studies of 
Widrow 2 and Press, Ryden and Spergel 3 have given much insight into the dynamics and 
evolution of finite thickness walls; however, it would be desirable to develop an analytic 
description of domain wall evolution which would preferably be model independent. In 
other words, we would like to obtain an approximate equation for a wall surface without 
having to calculate the detailed dynamics of the scalar field. 

We start by deriving the equations of motion for an infinitely thin domain wall (which 
turns out to be a generalization of the Nambu equation of motion for strings). Then we 
derive the leading order corrections to this equation for walls of finite thickness due to the 
extrinsic curvature of the wall and the wall’s self gravity. We do this by using an expansion 
scheme for the field theoretic equations of motion. For the non-gravitating case we check 
that the results obtained using the less satisfactory “effective action” method are the same. 
Next we comment on a misconception that has arisen in recent years regarding the rigidity 
of cosmic strings 4 . Returning to walls, we proceed by discussing the relative importance of 
the gravitational and field theoretic correction terms, and under which circumstances each 
are relevant. We then apply our equations to a few simple cases to illustrate the rigidifying 
(or otherwise) effect of the correction. We compare our results with those of Widrow, and 
discuss the breakdown of our approximation. 

A domain wall can occur during a phase transition into a broken symmetry state 
where the vacuum manifold has two or more disconnected regions. The wall becomes the 
boundary between regions lying in different vacua. To model the essential features of a 
thick wall, we will consider a scalar field with Lagrangian 


(we shall take the metric to be of the form (-, +, +, +)) and equation of motion 


v a v<v 


9V{4>) 

d<f> 


= 0 , 


2 


which contain symmetry. Vtf) is assumed to have two or more minima at 

non-zero values of Typically, V(*) is taken to be a Sine-Gordon or X<i> 4 potential. We 
will consider the tut where < * ± < * >o on either side of the wall, and therefore 

define the wall sur/nee, E, to be the three dimensional surface on which * = 0. In the 
typical example of a \<f>* wall in flat space, we have 


V{4>) = A {<f> 2 -v 2 ) 2 

4> 0 = T/tanh(\/2A7?z) 

Tab = -2A7 7 4 sech 4 (\/2A7 ? 2)( 5o6 - V 0 zV 6 z) . 


In this case, we see that 

_ 1 
v/2At7 

a = -\/2At7 3 
3 

represent the thickness and energy per unit area of this wall respectively; these will be our 
expansion parameters for the wall equations. More generally, if V represents the symmetry 
breaking scale, and A the self coupling of the field, then we expect the thickness of the wall 

always to be 0(1/ V^v) and its energy per unit area 0(V At? 3 ). 

This is the solution for a flat domain wall in a flat spacetime, in order to model a more 
general domain wall, we need to assume that the thickness of the wall is much smaller than 
both the extrinsic curvature of the wall and the horizon length (which is of order (Ga) 1 )• 
We may therefore split the field equations into their components orthogonal and tangential 
to E, the wall surface. This was done in ref. 5 using a Gauss Codazzi formalism, and so 
we wiU merely paraphrase the method here, and refer the reader to [5] for the details. 

We let n° be a unit geodesic normal vector field to E, and let z be the proper length 
along the integral curves of n°. Each z =const. surface then has unit normal n B , intrinsic 
metric h a b and extrinsic curvature K ab defined by 

hob = 9a b — n < I n b (2) 

Kab = h c a v c n b • 
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One can now use the Gauss Codazzi formalism to write the coupled Einstein-scalar equa- 
tions in terms of h a b, K ab and (f>: 

£nhab = 2 K a b 

(») (^ a ) 
C n K ah = Rab + 2K ac K£ - KK ab - 8TrGD a <t>D b <f> - 8tt GV(4>)h ab 

C n C n <f> + KC n <j) + D a D a <f> - = 0, (36) 

where D a and W R ab are the derivative operator and Ricci tensor of the wall surface. 

Now consider the new variables 


u = zjt X = </>/?7 

with e defined by (1), then C n — A ^ in this coordinate system. Moreover, since we can set 
A77 U(X ) = gC * <jU(X), where U(X) will now be 0(1), we may rewrite equations 

(3) as 


dfrq b 
du 


= 2eK ah 


dK ab 

du 


= -3irGcrU(X)h a b + e( W R ab + 2K ac K£ - KK ai ) - SttGctc 2 D a X D b X 


(4a) 


d 2 x idu(x) ax ~ 

du 2 2 dX + tK du + C DaD X ~ °* ( 46 ) 

We now write our solutions in terms of a power series in e. Clearly the zeroth order solution 
to equation (4b) is X = X 0 (u), where X 0 (u) is the flat space planar wall solution (tanh u 
for the \<f> 4 wall). We also see that h 0O * is constant, and K oab is given by 

K oah = tf o0 6(0) - 3nGcrh oa b J duU(X), 

which is just a constant term plus an odd function of u . 

Taking the difference of this last expression on either side of the wall yields the Israel 
equations for the exterior spacetime to the wall. On either side of the wall the Einstein 


4 



constraint equation W R + K ab K ab - K 2 = 0 is satisfied. Taking the difference of the 
constraint equation on either side of the wall and using the Israel equation we find that 
j£ o (0) = 0. There is another way to derive this result using only the scalar equation 
(4b). Multiply equation (4b) by dX/du and integrate from -oo to oo using the boundary 
conditions on X to obtain 

[°° du \{X') 2 K + tX'D a D a x\ = 0 . ( 5 ) 

J — oo 1 


Here a prime denotes derivative with respect to u. Now evaluating equation (5) to zeroth 
order in e we find 


/: 


iu(X'JK, = 0 
K 0 ( 0 ) = 0 . 


But this is just the wave equation for the wall surface. Thus, to zeroth order in wall 
thickness independent of the wall model, or any ansatz for the fields and even independent 
of gravity, the wall must obey a ^ambu 5 equation of motion. 

We will now investigate the case where the self gravity of the wall is negligible in the 
regimes of interest, in other words when {Git) <C e- Here, our equations of motion for 
X and h ab remain the same, but the equation for K ab simplifies. Since we are neglecting 
gravity, the metric is flat and one can show that ( *i2 ab = KK ab — K a c K cb • Now using 
equation (4a) and neglecting all terms proportional to G<r we find 

= —eK a c K c b . 

du 


The solution of this equation is 


K a b = K a b { 0 ) - euK a c {0)K c b 

Using the trace of this equation in equation (5) we find 


( 6 ) 


tf(0) r du{X') 2 + ef°° du f-u(X') 2 li: o 6(0)ii: b a + X'D a D a X 

J — oo V — OO 


= 0 


( 7 ) 
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We now find if(0) to second order in e by using an iterative procedure on equations (4b), 
(6) and (7). Evaluating equation (7) to zeroth order in e gives, as before, K o (0) = 0. It 
then follows from equation (6) that K 0 = 0. Now evaluating equation (4b) to first order 
and using the boundary conditions for X yields X x — 0. Equation (7) to first order then 
gives 1^(0) = 0. Now use equation (6) in equation (7) and evaluate to second order to 
obtain 

* ,(0) (/I i ' lX '°) + K ° tK °‘ K ‘‘ (/” iuu ' x °) = 0 

=> *,( 0 ) = -aK^KlKl, 


where 


a = 


f°° duu 2 X' 2 

J — OO 0 

f°° duX'1 

^ — OO 0 


> 0. 


( 8 ) 


Thus to second order in the wall thickness, the wall no longer obeys a Nambu equation, 
but acquires a term cubic in the curvature of the wall: 


K + ae 2 K° b K* c K' a = 0. (9) 

For the case of the \<f> 4 wall, one finds that a = ( 7r 2 — 6)/12 and for the Sine Gordon case 
a — 7t 2 /12. 

We now briefly remark upon walls with significant self-gravity. In this case, we have 
not only to worry about the behaviour of the wall surface itself, but also the way it affects 
the spacetime surrounding it. We therefore end up with two sets of equations which 
must be satisfied by the gravitating domain wall, the internal geometric equation, and the 
external Israel type equation. One can retain gravity in the above calculation if one wishes, 
and the geometric correction turns out to be similar and second order in e. However, the 
corrections to the Israel equations were worked out in [5] for the case of a A <f>* domain wall 
in vacuo, and were found to be first order in the wall thickness: 


[fc a fc] = —^TrGa-hab , (10) 

where 

cr = o (l + 7r eGcr (4 In 2 - j|)) , 
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and [kab] is the jump in the curvature across E, the surface representing the domain wall. 
For different models, the numerical factor in <r will change, but the form of the equation 
remains the same. Thus, for a strongly gravitating wall, the greater effect on the wall 
motion is due to its self gravity. We will return to this point later. 

Let us now consider the question of whether we can derive the equation of motion (10) 
in the absence of gravity, by an effective action argument. This would involve expanding 
the action 

s = J - V W] 

in terms of the wall thickness, and integrating out perpendicular to the wall. The method 
would involve using the coordinate system we have already set up, and proposing an ansatz 
for the scalar field: 4> — <f>o ( 2 ) "I" where <po satisfies 


and 6<f> satisfies 


a Vo _ av_ 

dz 2 d(po 


d 2 64> 

dz 2 


+ D a D a 6<p + K 


d<j>o 

dz 



0 . 


The action would then be given by 

s = + \J 

= J v=h{{- - n*o)][i + |z 2 (k ! - *D] + pwi} • (u) 

If one used the string argument, one would argue that 6<j> would have to be proportional to 
if, and thus that the final correction term would be a if 2 term. The equation of motion for 
this action is derived in the appendix. Noting that the action and therefore the equation 
of motion is valid only to order e 2 the equation of motion to order e 2 is 


D A D A X ja [l + |e 2 i?] + at 2 D B D A X oa D B D c X b 0 D A D c X ob = 0. 


where X 0 satisfies 


D a D a X 0 = 0. 
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Thus to second order the equation of motion is 

D a D a X, a = -ae 2 D B D A X oa D B D c X*D A DcX ob , (12) 

in agreement with equation (10). 

At this point we want to digress momentarily to address a misconception that has 
arisen concerning the rigidity of cosmic strings. The effective action for strings accurate 
to 0(e 2 ) was argued to be 9 


S = ~H J d 2 aV~h [l + e 2 A(K 2 + K\)] , (13) 

where A was some integral of the first order perturbations to the string fields. Using this 
action and the techniques of the appendix one can show that the equation of motion at 
zeroth order is just K 0 = 0 but this then implies that the first order perturbations to the 
fields vanish and hence A = 0. Thus the correct equation of motion for strings accurate to 
order e 2 is 

tfi(0) = K 2 (0) = 0. 

However, even if for some exotic type of string A was not equal to zero this would still 
be the correct equation of motion to second order. This is because for consistency of the 
expansion the zeroth order solution has to be used in the equation of motion to second 
order obtained by varying K 2 . 

Thus contrary to previous claims 9 to order e 2 strings are neither rigid nor antirigid. 
To determine their rigidity (or otherwise) the effective action to at least 0(e 4 ) must be 
obtained. 

Continuing with domain walls we now apply equation (10) to a few simple cases. 
Firstly we consider the collapse of a cylindrical wall. Taking cylindrical coordinates 
{t,p,<f>,z) for the 4-space and taking the radius p of the cylinder at time i to be P(t), 
the normal to the wall surface is n 0 = (1 - P 2 )~ 1 / 2 [V a p - PV a <]. Using equation (2) to 
find the extrinsic curvature of the wall we obtain 


K a b = (1 — P 2 )~ 3/2 


1 - P 2 
P 


^“•06 


Pu a u b , 
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where ip a and u a are given by i> a = a4> and u a = (1 - P 2 ) 1/2 [ V »< PV a p). Thus the 

equation of motion to zeroth and first order in e is 


r i 

= (1 _p2)3/2 + P(l-P2)!/ 2 


= 0 


i.e. 



This implies that 



(14) 


where Pi = — is the maximum radius of the cylinder and M is its mass per unit length. 

wviv- * S 2 tT<7 

This equation can be solved to obtain 


P — Pi cos 



where ti is the time at which P — 

To second order the equation of motion is 

K{ 0) = -SaKtKZKl 

Using the zeroth order equations of motion it is easy to show that the correction to the 
equation of motion is zero. Thus (14) is true even to second order. The equation of motion 

is valid provided 





as was previously suggested by Widrow 2 . Note that this does not agree with the naive 
expectation that the thin wall approximation is valid so long as the spatial radius of 
the cylinder is much greater than the wall’s thickness. The time components of K£ are 
important! 
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Now consider the collapse of a spherical wall of maximum radius Ri = ( — and 

mass M. In spherical coordinates (i,r,9,^) the wall’s radius is R(t) and the normal to 

the walls surface is n„ = (1 - - £v„i]. The extrinsic curvature of the wall is 

given by 


K a b = (1 - R 2 )~ 3/2 


R q b 


Ru a u b 


where q ab _ r 2 (V a #V{,0 -f sin 2 0V a <pV b <f>) and u a — (1 — B 2 ) _1 / 2 [V 0 < — _RV a r]. 
The equation of motion to zeroth and first order in e is 


K(0) = 


R o 


+ 


(1-R 2 ) 3 / 2 R 0 (l -Riy/> 


= 0 , 


which imphes that 

(1 -*>-(£)*■ 

This equation was integrated numerically using Romberg’s method of order 10 and the 
results is shown in figure 1. To second order in e the equation of motion is 


*2 + RUt) 


( 15 ) 


Note the correction term is positive and therefore finite thickness corrections lead to a 
slower collapse of the wall. 

Equation (15) can be integrated to obtain 


Rii = -y/l - x 4 e w( y -«-i) 


(16) 


where 


Figure 1 shows the result of numerically integrating (16) using a fourth order Runge-Kutta 
method for the case W=0.1. The equation (16) would be expected to be a poor description 
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of the collapse of the wall when 2 ^ 01 ? i e - 





) 


and invalid when \K£\ ~ e i.e. when y ~ W »/« [again in agreement with Widrow 2 ]. 
Note again that this is much sooner than naive expectations that would suggest that the 

approximation fails when the spatial radius of the sphere is O(e). 

In the case when ^ where R H ~ the horizon size associated with the 

wall, the most important corrections to the equation of motion arise due to self-gravity. 
Using (11) [see [5] for more details] the equation of motion for the collapse of a spherical 
wall in the vacuum can be written to first order in G<rRi as 


Riz 2 = 1 - 


(l-S(l-y 3 )) 2 


(17) 


where 

z = R 1 /R i , S = 2-K<rGRi. 

In figure 2 we show the results obtained by integrating (17) for the case 5 = 0.1. Gravity 
can clearly be seen to reduce the rate of collapse of the wall. 

Finally, it is easily checked that the exact solutions to the field theoretic equation of 
motion obtained by Vachaspati^ satisfy (10) to order e . 

To summarize, we have found the second order correction to the equation of motion 
of a domain wall, and exhibited that this gives a rigidity to the wall. We have shown 
that gravity, when important, has a stronger effect on the wall, also tending to make it 
more rigid. We estimated the breakdown of our approximation, which occurs typically well 
before the spatial radius of curvature becomes of the order of the wall thickness. Further, 
once this approximation does break down, corrections to all orders in thickness become 
important; in other words, only a field theory description will suffice to describe such a 
region. 

For the case of cosmic strings, we argued that there was no correction to second order 
in the width of the string, therefore strings are Nambu to a very high accuracy, unless they 
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exhibit regimes of high extrinsic curvature (such as cusps). At such places, only a field 
theoretic description will suffice. However, note that this breakdown will occur well before 
the spatial radius of curvature becomes of order the string width. 

Finally, Nielsen and Olesen were interested in a vortex solution exhibiting Nambu 
properties in order to get a phenomenological representation of a dual string. The lack 
of the second order correction term in the case of strings shows that the U yfj(l + j fc 2 )” 

lagrangian first suggested by Polyakov 8 , has no physical manifestation. 

Appendix. 

Here we will find the equations of motion associated with the action (11). In order to 
do this, we use an alternate definition of the curvature in terms of the coordinates of 2, 
which we will denote as £ A being the intrinsic coordinates of the surface 2: 

DaDbX* = Kab^* 1 (^ 1 ) 

We therefore need to find the equations of motion associated with the action 

s = -<rj d 3 ^[l + 0e 2 (D A D A X a ) 2 - A e 2 D A D B X a D A D B X a }. (A2) 

where A = a/2. In varying the action with respect to X\ we must remember that both 
the metric, h AB , and the connection, h E B depend on X a . For the metric we have: 

Sh AB = 6X° A X a>B + X° A 6X a>B (A3) 

Therefore, varying 5 with respect to X a gives three pieces The first term is the variation 
of the volume element, this gives, upon integration by parts, 

6Sl = a J d 3 ty/=h6X a D A [D A X a [l+0e 2 (D c D c X b ) 2 - A e 2 (D B D c X a ) 2 }] (A4) 

The other two terms contain the variation 

6[h AC D A D B X a ) = -h AE h CD 6h ED D A D B X a + D c D B SX a - h AC 6h E B X* E (A5) 
the last term of which vanishes when contracted with DDX a , which is normal to the wall. 
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Substituting from (A3) and (A5) and integrating by parts as necessary gives 
6S 2 = cr/3t 2 J d 3 (^h^D c [D B D B X b D A D c X b D A X a ] + 2(D A D A ) 2 X a }6X a , ( A6 ) 
and 

<S, = -<rAe 2 [dHV^h{iDc [D c D B X l D A D B X t D A X a ] +2 D b D a D a D b X.}6X°. 

(At) 

Now, the Gauss Codazzi equations imply 

R ac = KK ac - K ab K b = D B D B X b D A DcX b - D A D B X b D B D c X b (AS) 


(A9) 


(A10) 


and hence 

6S 3 = -<rAe 2 J d 3 is/-h{AD c [D B D B X b D A D c X b D A X a ] - 2 D B [R AB D A X a ] 

+ 2D A D A D B D B X a }6X a . 

Thus, imposing ^ = 0, (A4), (A6) and (A9) give as the equations of motion for the wall: 
D a [D A X a [l+(3t 2 {D c D c X b ) 2 - A t 2 (D B D c X a ) 2 }] +2(0 - A )e 2 (D A D A ) 2 X B 
+ 4({3-A)t 2 D c [D B D B X b D A D c X b D A X a ] -2A e 2 D B [. R AB D A X a ] = 0. 

Now, from our previous work, we know that the first order corrections vanish, and hence 
0 = A. This can also be deduced if we use the zeroth order form of the above equation 
( D A D A X a = 0) and input this information back into the argument at the stage of ‘solving’ 
for the field perturbations. (One may regard this as slightly circular.) However, using this 
simplification one obtains 

D A D A X a [ 1 + Ac 2 *] + 2A e 2 D B D A X a D B D c X b D A D c X b = 0. (All) 

Which is the equation of motion we wished to derive. 
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Figure Captions 

(l)The collapse of a spherically symmetric domain wall to zeroth order (solid line) 
and second order (long dashed line). The horizontal dashed line corresponds to y = 



(2)The collapse of a spherically symmetric domain wall for the case S=0.1 to zeroth 
order (solid line) and second order (dashed line) 

References 

[1] C.Hill, J.Fry, D. Schramm, Phys. Rev. D39 3571 (1989). 

A.Stebbins, M. Turner, FERMILAB-PUB-88/170-A 

[2] L.Widrow, Phys. Rev. D39 3571 (1989). and Phys. Rev. D40 1002 (1989). 

[3] W. Press, B.Ryden, D.Spergel, NSF-ITP-89-51 

[4] This has independently been noted by J. Simon, UCSB-89-50 

[5] D.Garfinkle, R.Gregory, FERMILAB-Pub-89/255-A 

[6] K.Maeda, N.Turok, Phys. Lett. 3B 376 (1988). 

R.Gregory, Phys. Lett. 199B 206 (1988). 

[7] T.Vachaspati, Vachaspati, TUTP 89-7 

[8] A. Polyakov, Nucl. Phys. B268 406 (1986). 


14 



initial 




initial 



